Scaling and exact solutions for the flux creep problem in a slab superconductor 
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The flux creep problem for a superconductor slab placed in a constant or time-dependent magnetic 
field is considered. Logarithmic dependence of the activation energy on the current density is 
assumed, U — Uo ln( J/ Jc), with a field dependent Jc. The density B of the magnetic flux penetrating 
into the superconductor, is shown to obey a scaling law, i.e., the profiles B{x) at different times t 
can be scaled to a function of a single variable x/t'^. We found exact solution for the scaling function 
in some specific cases, and an approximate solution for a general case. The scaling also holds for a 
slab carrying transport current I resulting in a voltage V (x P , where p ~ 1. When the flux fronts 
moving from two sides of the slab collapse at the center, the scaling is broken and V(I) crosses over 
to 1/ oc /^o/^T. 
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I. INTRODUCTION 



Thermally-activated hopping of flux lines between pin- 
ning sites, or flux creep, controls both the magnetic and 
transport properties of superconductors under various 
external conditions. In particular, it is responsible for 
the frequently observed fast magnetic relaxation and the 
highly nonlinear local voltage-current curves. Flux creep 
is specifically pronounced in high-temperature supercon- 
ductors because there the pinning energijea. are small, 
while the operation temperatures are high.ErB 

In the majority-of theoretical and experimental stud- 
ies of flux creeptTLl a superconductor is flrst placed in an 
increasing or decreasing field, thus acquiring a non-zero 
magnetization, M . The field is then kept constant and 
the relaxation of M with time is examined. In the present 
work we address a different problem - flux penetration 
into a non-magnetized superconductor, M{t = 0) = 0. 
After turning on a constant or time-dependent magnetic 
field, a flux front propagates from the surface with some 
time-dependent velocity. In the case of a long slab in par- 
allel field two planar flux fronts propagate from each side. 
Absence of a characteristic spatial scale suggests a pos- 
sibility for scaling solutions, i.e., the flux density proflle 
B{x) at different times t is a function of a single variable 
x/t^ , where /3 is a constant. Indeed, such scaling was 
demonstrated in Ref. |^ for an instantaneous turn-on of 
a constant applied field, and assuming a superconductor 
characterized by a logarithmic dependence of the pinning 
energy on the current density. 

In this work we have sought for the whole class of flux 
creep problems having solutions with a single-parameter 
scaling. This class turned out to be a lot broader than 
considered in Ref. ^, insofar as it extends to situations 
with time-dependent applied fields, Ba{t) cx , and a 
general field-dependent critical current density Jc{,B). 
We also report scaling solutions for the flux creep in a 
superconductor carrying a transport current. There the 
flux creep manifests itself in a number of experimental 
observations like a relaxation in the resistance, and that 



the voltage acij£iss-.the sample depends on the current 
sweep rate, etcljHj 

The paper is organized as follows. In Sec. II the ba- 
sic equations are formulated. Section III brings out the 
scaling properties of their solution. In Sec. IV two cases 
allowing an exact analytical solution are considered. The 
subsequent Sections report the implications of our results 
for measurable quantities - J, E distributions, mag- 
netization and voltage. Finally, the main conclusions are 
summarized. 




FIG. 1. Superconductor slab in an applied magnetic fleld 



II. FORMULATION OF THE PROBLEM 

Consider a semi-infinite superconductor (filling the 
half-space a; > 0) placed in an external magnetic field 
Ba{t) directed along the 2:-axis, see Fig. |l|. The position- 
and time-dependent fiux density (or magnetic induction) 
and shielding current in the sample are denoted as -B(a;, t) 
and J(x,t), respectively. Both J and the electric field E 
are directed along the y axis. 

From the Ampere and Faraday laws one then has, 



^oJ = —dB/dx , 
dB/dt= -dE/dx. 



(1) 
(2) 
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To describe the superconductor we assume that E = 
vB, where v = fo exp(— JJ/ZcT) is the velocity of the 
thermaUy-activated vortex motion over the barrier U. 
With a logarithmic dependence U{J) = Uo\n{Jc/J) it 
follows that 



E{J,B)=vo\B\ |J/J,|"sgnJ, 



(3) 



where the exponent n = Uo/ksT depends on the temper- 
ature T. The critical current density Jc usually depends 
on \B\, and in the bulk part of|-the paper it will be as- 
sumed that this is described by,tHl 



Jc{B) = J,o \Bo/B\ 



(4) 



To avoid parameters with inconvenient dimensions wc in- 
troduced here two constants, Jco and Bq, a characteristic 
currentpdensity and induction, in addition to the expo- 
nent 7.t3 With 7 = one has a constant Jc (the Bean 
model) , while for 7 = — 1/n one has an E — J relation, 



E = vqBq \ J/Jco\" sgn J, 



(5) 



which is B- independent. Such an E(J) has been used by 
many authors since it allows a significant simplification 
of the analysis. It is justified when the magnetic field 
inside the superconductor is essentially uniform, in par- 
ticular, when a small perturbation 6B{x, t) is added to a 
superconductor cooled in a large dc field. 

The electrodynamic problem defined by Eqs. (|l|)-(Q) 
can be reformulated as a non-linear diffusion equation 
for the flux density. 



dB 
'dt 



{(J-oJcqB^ 
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Similar equations can be written also for J and E. 

We will first consider a totally flux free superconductor 
which at i > experiences an increasing applied mag- 
netic fleld given by 

^a(t) = So(t/Tr , a>0, (7) 

where t is another constant. For a = and a = 1 this 
describes an instant field step and a linear ramping up, 
respectively. By introducing the dimensionless variables 

fJ-oJcO r t , J E 



Bo 



~ t J 

T JcO 



VqBq 



(8) 



and removing the redundant parameter definition by 
choosing 



the Eq. 



Bo = fioJcQVoT , 

acquires the formlli 



dt di V ^ 



db 
dx 



The boundary condition becomes 
at the superconductor surface. 



(9) 



(10) 



(11) 



III. SCALING 



The flux creep problem Eqs. ( p^ and (^Tj) is solved by 
writing the flux density in the scaling form 



5(£,f) -t"/(0, ^ = ii- 



with 



1 + an{l + 7) 
1 + n 



(12) 



(13) 



By substitution one flnds that the scaling function /(^) 
satisfies the differential equation 



(14) 



where it was used that / > 0, and /' < 0, i.e., the flux 
density decreases monotonously as one moves away from 
the surface. The boundary conditions become 



/(0) = 1, /(eo) = 0. 



(15) 



where the last one arises from the physical requirement 
that the scaling function vanishes at the flux front, which 
is located at ^ — £,o- From Eq. (12) it follows immediately 
that the flux front advances with time according to 



Xo 



(16) 



For the case a = (constant applied field) one has 
/3 = l/(n-|-l), while a larger a increases the exponent (3. 

Using Eqs. (|l|) and (^) one finds that also the current 
density and electric field have scaling properties, namely 

e{i,i)^i"+P-' feiO, /. = |/r/'+^"- (17) 

The scaling behavior of the ac losses is found from a sim- 
ilar analysis of the product j(x, t)e{x, t). 

When the magnetic field is instantly turned on and 
then kept constant, (a = 0), exact scaling relations are 
obtained even when the power-law Eq. (Q) is relaxed and 
replaced by any Jc{B). Considering here 



Jc{B) - 



s{B/Bo) ' 



(18) 



where s is a general function, a diffusion equation for the 
flux density similar to Eq. ( p^ can be derived giving. 



db 
'dt 



d_ 

dx 



\b\s-(b) 



db 
dx 



db 
dx 



(19) 
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where the same dimensionless variables Eqs. (||) and (||) 
are used. The scahng form of the flux density is 
still applicable where now b(x,i) = The 
equation for the scaling function /(^) is changed to 



e/7(n +!) = [/ i/'n' 



(20) 



while the boundary conditions ( |l5|) are the same. The 
flux front, current density and electric field are now given 
by 

- h ^ l/r / , (21) 

hence, also they possess scaling properties in this case. 

/^ocB"" y^ = const y^ocB T General y^(B) Jci^) 







Sjj=const 
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FIG. 2. Summary of our results. The diagram illustrates 
where the flux creep problem has scaling solutions. Full ana- 
lytical solution is found for two particular cases. 

As we have formulated the flux creep problem, the wide 
variety of possible situations is represented by the set of 
three independent parameters, n, 7 and a. For a given 
parameter set the scaling function /(^) can be obtained 
by solving the Eqs. ( p^ ) and ( |20| ) numerically. Interest- 
ingly, one may in two special cases find exact analytical 
solutions, as shown in the next Section. Shown in Fig. ^ is 
a summary of our findings regarding scaling solutions for 
the flux creep problem in various situations. From this 
"Jc{B) — Ba{t)" diagram it is seen that scaling holds for a 
majority of conceivable conditions. Note that even if the 
scaling function itself is unknown one may test predicted 
scaling behaviors experimentally, e. g., by analyzing flux 
density profiles B{x, t) measured at different times. Plot- 
ting B{x, t)l Ba{t) versus x/xo(t) is expected to result in 
a collapse into a scaling function /(^). Other practical 
consequences of the scaling concerning the magnetization 
and voltage behavior are discussed below. 



IV. EXACT SOLUTION FOR PARTICULAR 
CASES 

A. Constant field, B-independent E(J) 

The solution for this case will be worked out in con- 
siderable detail in order to illustrate some basic features 



of the scaling properties of the physical quantities. Sup- 
pose that at i = a constant applied magnetic field, 
-Bo, is turned on, i. e., one has a = 0. Furthermore, 
let 7 = —^/n, which leads to the commonly assumed 
current- voltage law, Eq. (||). From Eqs. (|l2[) and ( p^ ) 
the scaling law has then the form 

b{i,i) = f{0, £.^xi-^/(^^+^\ (22) 

and Eq. (|lj) for the scaling function reduces to 

[^-n(n + l)|/'rV"] /'-O. (23) 

This equation has two real solutions for /' consistent with 
dB/dx < 0, namely 



-/' = 



2n 



7 TTT (^0 

(n + 1) 



i/("-i) 







(24) 



The parameter ^0 is an integration constant which is de- 
termined from the following interpretation of Eq. (|24|). 
Noting that since —f'{£,) directly represents the spatial 
profile of the current J{x,t), the solution with /' = 
is identified as the one valid for the Meissner-statc re- 
gion where also 6 = 0. The point ^ = ^o, then defining 
the flux penetration front, is determined by /(^o) = 
using the solution valid for the penetrated (mixed-state) 
region. We find 
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n 



1 



F(l) 



l-n 



l/(n+l) 



1 

where the function F is defined as 

^(l-y2)i/("-i'dy 





F{z) 



(25) 



(26) 



The value F{1) can be expressed through F- functions as 

'3 1 



F(l) 



r 



/r 



2 V" - 1/ \2 " - 1 
The result for the fiux density profile becomes 

F{i/ioii)) 



b{x,t) ^ 1 - 



Fil) 



< X < Xq , 



(27) 



(28) 



where the position of the flux front, xq, is moving with 
time according to 



Xq 



(i)=eoP/(n+l) . 



(29) 



It then also follows that the normalized current den- 
sity j = —db/dx, and electric field in the flux penetrated 
region, < x < xq, are given by 



1 



1 



e{x,i) = j{x,i)" . 



Xo (t) 



(30) 



3 



The differential resistivity dE/dJ oc therefore 
varies in space as a parabola having a maximum value 
at the edge. 

Figure |3| shows the distribution of the magnetic induc- 
tion, current density and electric field for n — i and 11. 
As the time increases, the flux penetrates deeper into 
the sample accompanied by a smaller slope of the flux 
density profile, i. e., a decay of current density. The elec- 
tric field profile follows the behavior of J(x), although 



decreasing much more rapidly due to the strongly non- 
linear E{J) law. As the power n increases the field and 
current distributions are seen to become more linear. In 
the limit n = oo one has F{z) = z and = Ij f-nd hence 
io(t) — 1. The temporal dependences then vanish, and 
the behavior reduces to 6 = 1 — a;/(5, = 1, and e = for 
< X < where 8 — B^j [IqJcq. Thus, the well-known 
results of the Bean model are reproduced. 
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FIG. 3. Distributions of the normalized flux density, current density and electric field after a sudden turn-on of a constant 
applied magnetic field. A B-independent E{J), i.e., 7 = — 1/n, is assumed, and the plots are made for n = 3 and 11. The 
arrows indicate time direction, and the curves correspond to t =0.05, 0.2, 1, and 3. 



B. Linearly increasing field, constant Jc 

For the linear sweep case, Ba{t) — Bat {Ba = const), 
and a constant Jc, i. e., a = 1 and 7 = 0, the flux creep 
problem has a surprisingly simple solution. It follows 
from Eq. @ that 

B{x,t) ^ Bat[l~x/xoit)], x<xo, (31) 

and the region x > xq is flux free. The flux front propa- 
gates linearly in time, 

D / T \ + 



Remarkably, this solution coincides with the Bean model 
behavior where in the penetrated region the flux profile 
is a straight line, and the current density is uniform, 

J{x) ^ Jc iBa/fioJcVoj = const. (33) 

The electric field E{x) becomes proportional to B{x), 

/ ■ \ n/{n+l) 

E{X, t) = Vo [Ba/tJLQJcVQ) B{x, t) , (34) 

hence, also this quantity is linear in both x and t. All 
these profiles are presented in Fig. ^(b). For larger field 
sweep rates, Ba, the current density and electric field 
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are higher, while the penetration depth at a given Ba is 
smaller. 

The results of this subsection, unlike the results of the 
rest of the paper, can be used also to describe the flux dis- 
tribution at large enough fields when the flux penetrates 
the whole slab. The current then flows in the opposite 
directions in the left and right halves of the slab, but with 
the same density as defined by Eq. (|3^). The magnetic 
field in the left half is given by a linear profile defined 
by Eq. (^, where xq is still given by Eq. (|3^), though 
it does not have the meaning of the flux front position 
anymore. In the right half B{x) is a mirror image of that 
in the left half. 



V. SUMMARY OF B, J, AND E 

The evolution of the flux density, current density, and 
electric field distributions found numerically is presented 
in several graphs. Figure ^ shows results for the applied 
field linearly increasing in time, while Fig. ^ as well as 
Fig. 1^ - for an instant turn-on of a constant field. From 
these graphs the following conclusions can be drawn. 

(i) The flux density profile is found convex for 7 > 0, 
i. e., for constant or monotonously decreasing with 
B. It also means that dJ/dx > 0, i. e. the current 
density increases as x approaches the flux front. For 



_B-independent E{J) the picture is the opposite; |J| is 
maximum at the edge, and the B(x) profile is concave, 
see Fig. || and Fig. ||(a). This behavior is expected for a 
field-cooled superconductor when a small additional field 
is applied. When comparing these results to experimen- 
tal data one should keep in mind that the shape of the ob- 
served B{x) profile can easily be affected by demagneti- 
zation effects. However, the sign of dJ/dx for the current 
density profile seems to be robust, which is confirmed, bjt 
the flux creep simulations for a thin film geometry.EjEEl 
Schematic J{x) profiles for the most important cases are 
shown in Tab. 1. 

(ii) Even when Ba increases linearly with time, the flux 
penetration can proceed with acceleration or deceleration 
depending on the Jc{B) dependence. This is illustrated 
in Fig. ^ where the panels (a)-(d) are ordered according 
to the Jc{B) behavior, and the time intervals between 
the curves are equal. Only in (b) is the flux front mov- 
ing with a constant velocity. In (a) the speed is slow- 
ing down, while in (c) and (d) it is increasing. In the 
latter case, ramping of the applied field leads to an ef- 
fective reduction of Jc, and hence to larger electric fields 
and faster penetration, (a), (b), and (c,d) correspond to 
P < 1, (3 — 1 and /3 > 1, respectively, where the expo- 
nent j3 given by Eq. (p^) c ontrols the growth of the flux 
penetrated region, Eq. (|16D. 



TABLE I. The shape of the current density profiles and the parameter Co characterizing the flux front dynamics, Eq. (p^, 
in different regimes. The parameter a describes the approximate solution Eq. for the scaling function /. 



n 


B- independent E{j) 
(7 = -1/n) 


jc — const 
(7 = 0) 


j4B) <x 1/B 
(7 = 1) 


The Kim model 
Jc = Jco{l + B/kBo)~^ 




Constant field, Ba = BoO{t) {a = 0) 




J 


^ X 


J 




X 


J 




1. 


J 




X 




see Fig. |^ 


see Fig. p|(a) 


see Fig. ^(b) 


see Fig. p|(c,d) 


3 


exact solution, 
Eq. @ 


Co = 1.521, a = 0.105 


Co = 0.881, a = 0.032 


Co = 2.044, a = 0.42 for fc=l, n=3 
Co = 1.702, a = 0.25 for fc=3, n=3 


7 


Co = 1.346, a = 0.045 


Co = 0.724, a = 0.014 


15 


Co = 1.212, a = 0.021 


Co = 0.628, a = 0.007 


33 


Co = 1.119, a = 0.010 


Co = 0.568, a = 0.004 




Linear increase, Ba = Ba t {a — 1) 




J 




J 






J 




X 


NO 
SCALING 


1 X 




X 


see Fig. ^(a) 


see Fig. ^(b) 


see Fig. ^(d) 


3 


Co = 1.468, a = 0.016 


exact solution, 
Eqs. (|l[||) 


Co = 0.514, a = -0.004 


7 


Co = 1.164, a = 


Co = 0.503, a = -0.002 


15 


Co = 1.071, a = 


Co = 0.501, a = 0.000 


33 


Co = 1.031, a = 


Co = 0.499, a = 0.001 
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FIG. 4. Distributions of the normalized flux density, current density and electric field at equidistant times, t =0.5, 1, 1.5, 
and 2. The applied field is linearly increasing, and n = 3. (a)-( d) diff er in Jc{B) dependence: (a) B-independent E{j), Eq. (|^; 
(b) constant Jc {Jc ~ 2Ba/fJ^oVo), a case solved exactly in Sec. [VB; (c) and (d) Jc decreases with increase of \B\. The arrows 
indicate the time direction. 
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FIG. 5. Distributions of the normalized fiux density, current density and electric field at different times, t =0.05, 0.2, f , and 
3. The applied field is turned on at f = and kept constant, and n = 3. (a)-(d) differ in Jc{B) dependence: (a) constant Jc; 
(b) Jc = Jco/b; (c) and (d) the Kim model, Jc = Jco/(l + 6) and Jc = Jco/{l + b/3), respectively. The arrows indicate the time 
direction. 
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(iii) The electric field has always a maximum at the 
edge and decreases monotonously with which follows 
from dE/dx = ~B < 0. Physically, this occurs simply 
because all the vortices enter the slab through the edge 
region where the flux motion is thus most intense. The 
second space derivative of E is proportional to dJ/dt, 
and one can find different curvatures of the £'-field pro- 
file for different cases. 



VI. MAGNETIZATION 

Since the nonlinear flux diffusion has a sharp front lo- 
cated a flnite distance from the surface, the solution ap- 
plies also to the case of a slab of finite width 2w provided 
w is larger than the penetration depth, xq ■ The flux then 
penetrates the slab from both sides with non-overlapping 
profiles. 

The time-dependent magnetization of the slab is given 

by 



B(x)dx , 



(35) 



where B is non-zero only in the region, < x < XQ{t). 
Substituting here the scaling law, Eq. (p2), one obtains 



fioM{t) = Bait) 



I- A 



Xo{t) 



0<xo<w, (36) 



where A is a number equal to the average value of the 
scaling function in the penetrated region: 



A 



1 



(37) 







For the exact solution considered in Sec. IV A one finds 
A = (n- l)/[2ni^(l)] < 1/2. 



For the exact case in 
Sec. p^VB| where the B-profile is linear, one has A = 1/2. 



For fiux density profiles with a convex shape, which are 
found when Jc decreases with B, one has A > 1/2. 

Expressing the time-dependence explicitly, the magne- 
tization can be written as 



fioM{i) ^ Ba{i) [l-A{i/i*y 



<t <t* 



(38) 



where t* is the time it takes for the fiux to completely 
penetrate the slab, i. e. i* — (w/S^o)^. For the fully- 
penetrated state, i > i*, the time-dependence of M is 
essentially different, and one expects a kink in the mag- 
netization relaxation rate, as first predicted in Ref. W 



VII. TRANSPORT CURRENT AND 
NON-STATIONARY V(I) 

Most results found in the paper are also relevant to a 
superconductor carrying a transport current. The flux 
penetration into a slab biased with transport current is 



governed by the same equations as those for a slab placed 
in an applied magnetic field. The difference is that for 
the transport current case the current fiows in the same 
direction on both edges of the slab. However, this is 
not important as long as the two penetrating flux fronts 
do not meet in the slab center. The boundary condi- 
tion B{x = 0) = Ba for the applied field case should 
be replaced now by B{x = 0) = /xo//2, where / is the 
transport current per unit height of the slab. 

It is interesting to analyze the voltage V measured on 
the superconductor carrying a transport current /. It is 
obvious that V is not a unique function of / but also 
strongly depends on how / changes with time. In a con- 
ventional experimental setup the voltage contacts are at- 
tached to the superconductor surface, thus the measured 
voltage is determined by the electric field in the surface 
layer only, V{t) = E{x ~ 0, t)L, where L is the distance 
along y between the contacts. 

The electric field at the edge has a power-dependence 
on time given by Eq. ([T^), since /e(0) is just a number. 
Thus, from Eq. ( p^ one obtains the following power-law 
voltage-current relation, 



VocIP, p=l 



1-1 
a 



(39) 



Note that the exponent p of the integral V{I) can be 
much different from the exponent n characterizing the 
local E{j). Actually, p is only weakly sensitive to n, espe- 
cially for large n. Moreover, increase in n can sometimes 
lead to a decrease in p. 

For large transport currents when the flux penetrates 
the whole sample the present analysis is not valid. Then 
the current density will be distributed over the slab more 
or less uniformly, and one can expect that the inte- 
gral voltage-current curve will reflect the local one, i. e., 
V oc J""*"^. The crossover between low-current and high- 
current parts in V(I) has been reproduced by numerical 
simulations in Ref. H. This crossover should be accessible 
experimentally since the exponent p for an incomplete 
flux penetration is of the order of unity, while the expo- 
nent n for the full penetration is temperature-dependent 
and can be very large. 

The exponent a in Eq. ( |39| ) dcflnes the time depen- 
dence of the transport current, I(t) cx t". For larger a 
V{I) becomes steeper until the exponent p saturates. For 
a = 0, i. e. when the current is turned on and kept con- 
stant, the voltage decays with time. This voltagMclax- 
ation has been earlier observed expecipjentallyoilil and 
reproduced by numerical simulations Jj'lj It is accompa- 
nied by relaxation of current distribution, as illustrated 
by Fig. ^. According to Eq. (|T7|), the voltage decay is 
described by oc t""/'""*"^', in particular y cx 1/t at low 
temperatures when n is large. 

Remarkably, for small enough a the voltage will de- 
crease with increasing current. For example, for a = 1/3 
and 7 = 0, one obtains 
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V cx /■ 



(40) 



in particular, F cx 1// at low temperatures. Physically, 
a decrease of voltage is related to the same relaxation 
process which takes place for any external conditions. 
Increase of / with time supplies more current to the su- 
perconductor and tends to increase the voltage. If / in- 
creases slow enough, then the first process is dominant, 
and the measured V{I) curve should give voltage decreas- 
ing as the current increases. The evolution of flux, cur- 
rent and electric field distributions for this interesting 
case is shown in Fig. ^. 



— 0.5- 




FIG. 6. Distributions of the normalized flux density, cur- 
rent density and electric fleld at different times, t =0.5, 1, 2, 
and 4. The transport current increases with time as 7 oc t^^^, 
while Jc = const and n — 3. The voltage, which is pro- 
portional to e{x = 0), decreases with time, as shown by the 
arrows indicating the time direction. 

For the experimentally most relevant case of a linearly 
increasing current, I ~ It, (i.e. a ~ 1) one obtains from 
Eqs. and (|) 



V 



Mo^^oi|/'(0)|" 



1^+^. (41) 



2 \^2i+7«oJcoBj; 
For a constant Jc, or 7 = 0, we come to the exactly 



solvable case considered in Sec. [VB, and the expression 
further simplifies to. 



V = 



I 



n/(n+l) 



2JcVo 



I. 



(42) 



In this case the superconductor behaves like an Ohmic 
conductor. Deviations from the Ohmic behavior can be 
caused by a i?-dependence of Jc- In particular, if Jc de- 
creases with B, the exponent of the V{I) curve becomes 
larger than unity, see Eq. (|4^). Meanwhile, at small / 
the superconductor always behaves Ohmically since the 
self- field is small and the Jc{B) dependence can be ig- 
nored. These results are in agreement with numerical 
simulations reported in Ref. 

We also note that when the transport current is 
ramped faster, the voltage at a given current is larger, 
V (X i "/("+!). This observation is in agreement with 
results of niuperical simulations and experiment on Bi- 
based tapes. □ 



VIII. APPROXIMATE SOLUTION 



In a general case Eqs. (|1J) and ( |20| ) with boundary con- 
ditions (^5|) cannot be solved analytically. Near the flux 
front / ^ 0, and Eq. (14) reduces to (const -I- /3Co)^^" — 



P |/'|. Therefore, the behavior of the scaling function / 
is determined only by the Jc{B) law at small B, and 



/(e) « (eo - 



1/(1+7) 



(43) 



This result holds true for any Jc{B) which has asymp- 
totic behavior Jc oc B^'' at i? 0, e.g., 7 = for the 
Kim model. 

Surprisingly, a very good approximate solution for / 
in the whole region < ^ !i Co is given by the expression 



/(e) = (l-e/eo)'/^'+"'(l + <) 



(44) 



Values of ^0 and a for several common dependences 
Jc{B), including the Kim model, and values of n have 
been found numerically and listed in Tab. 1. Expressions 
for ^0 E^nd a in a general case can be found by substituting 
/ from Eq. (^) into Eq. (|lj) and analyzing expansion in 
powers of (^0 ~ C)- Then, one obtains: 



1 2n(l +p)(p + an) 
a 



p'^{p — a) 



e, 



n 



1 



1, p= 1/(1 + 7), 
(l + a)"(i+7) 







l + an(l+7) (1 + 7)^ 



(45) 
(46) 



Since the scaling function /(^) usually has a very simple 
shape, we find that the approximate expressions fit the 
exact solutions with a good accuracy, e. g., the deviation 
is less than 1% for all cases shown in Figs. ^ and ^. 
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IX. CONCLUSIONS 
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The propagation of magnetic flux into a slab supercon- 
ductor has been considered using the flux creep approach 
with a logarithmic current dependence of the activation 
energy. The dynamic behavior was found to possess scal- 
ing for a general Jc{B) when a constant magnetic field is 
suddenly applied to the superconductor, and for a power- 
law Jc{B) in case of an applied field ramped up with a 
general power-dependence on time. For two particular 
cases of the creep problem an exact analytical solution 
could be found. 

The main results obtained in this work are: 

1. The flux density profiles at different times follow 
the scaling law, B{x,t) — Ba{t)f{xt~^). Similar 
scaling applies to the current density and electric 
field profiles. 

2. The flux density profile is convex for penetration 
into zero-field-cooled slab, and concave for a slab 
cooled in a large field. 

3. At constant Jc and linearly increasing Ba the Bix) 
and J{x) profiles at any time coincide with the 
Bean-model profiles. 

4. The flux front position is a power function of time 
given by Eqs. (jl^) and (p^). The front moves 
through the slab with an increasing or decreasing 
velocity depending on the material's Jc{B). 

5. The explicit time dependence of the magnetization 
is found. 

6. For a partially penetrated slab carrying a transport 
current /, the voltage is a power of /, with an 
exponent different from that of the local E{J) re- 
lation. A pronounced crossover in the V{I) curve 
at the point of full penetration is predicted. 

7. For a small transport current increasing linearly 
with time, the Ohmic behavior cx /, is found. 

8. For a stationary transport current the voltage de- 
cays as F ~ 1/t. 

9. An increase of transport current can be accompa- 
nied by a decrease of voltage, in particular, V ^ l/I 
when / cx t^^^. 

All the conclusions can be tested experimentally: (1)- 
(4) by spatially-resolved techniques, and (5)- (9) by inte- 
gral measurements. Our results presented by Eqs. (|3^ ) 
and ( |39| ) allow a new method to infer the material prop- 
erties such as local E{J) or Jc{B) characteristics from 
integral measurements of magnetization and voltage. 
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